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SECTION  I 
INTRODUCTION 

Aircraft  with  a strategic  mission  in  a general  war  or  in  any  conflict 
involving  nuclear  weapons  could  be  subjected  to  various  environments  generated 
by  the  detonations  of  nuclear  weapons.  During  a representative  mission,  the 
aircraft  could  be  targeted  by  eneiry  submarine-launched  ballistic  missiles 
(SLBMs)  during  base  escape,  it  could  penetrate  late-time  radioactive  clouds 
generated  by  detonations  of  enemy  nuclear  weapons  at  missile  sites  or  other 
hard  targets,  it  could  be  targeted  by  surface  to  air  missiles  (SAMs)  and/or 
by  air  to  air  missiles  (AAMs),  it  could  be  exposed  to  environments  generated 
by  collateral  weapon  detonations,  and  it  could  be  subjected  to  ionizing 
radiation  during  damage  assessment  surveys. 

At  the  present  time  there  are  analytical  tools  available  and  documented  to 
aid  the  analyst  in  investigating  the  majority  of  the  postulated  situations. 
Nuclear  blast  and  thermal  environments  can  be  predicted  with  reasonable 
accuracy  using  computer  codes  such  as  SPARK*  (ref.  1)  ano  SNARE*  (ref.  2).  The 
aircraft  response  can  be  estimated  using  VIBRA-4  (ref.  3),  a nuclear  gust 
response  code;  NOVA*  (ref.  1),  a nuclear  overpressure  response  code,  and  TRAP* 
(ref.  2),  a nuclear  thermal  response  code.  Electromagnetic  pulse  (EMP)  environ- 
ments can  be  estimated  using  the  EMP  Phenomenology  Handbook  (reference  4) 
or  EMP  environmental  prediction  codes.  The  system  EMP  response  is  not  straight- 
forward, but  efforts  are  underway  at  the  Air  Force  Weapons  Laboratory  to  develop 
the  technology  base  required  for  such  analyses.  Reference  5 can  be  used  to 
estimate  system  responses  to  EMP. 

Prompt  nuclear  radiation  environments  can  be  predicted  using  SMAUG  (ref.  6) 
and  initial  nuclear  radiation  environments  by  using  FIREFLY**  (ref.  7). 


*SPARK  and  SNARE  are  environmental  subroutines  contained  in  the  NOVA  and 
TRAP  computer  codes. 

**Prompt  nuclear  radiation  environments  are  those  generated  by  detonation  and 
are  associated  with  pulse  durations  of  a few  microseconds  in  duration.  Ini- 
tial nuclear  radiation  environments  include  fission  product  gamma  radiation 
and  are  associated  with  times  up  to  one  minute  after  detonation.  Initial 
nuclear  radiation  environments  include  the  prompt  environments. 
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Ionizing  dose  accumulated  during  radioactive  cloud  penetrations  can  be  esti- 
mated usino  the  results  presented  in  reference  8.  Response  to  nuclear 
radiation  generally  is  associated  with  electronic  equipment  and/or  the  aircrew 
and  can  be  estimated  using  available  techniques.  For  example,  equipment 
response  can  be  estimated  using  manual  circuit  analysis  techniques  or  computer 
codes  such  as  SCEPTRE  (refs.  9 & 10)  and  human  response  can  be  estimated  using 
techniques  and  information  from  such  publications  as  references  11,  12,  and  13. 

However,  there  is  one  situation  in  which  the  tools  noted  above  are  inade- 
quate, namely  damage  assessment  surveys.  A part  of  the  mission  of  some  air- 
craft, e.g.,  manned  bombers,  is  the  surveying  of  targets  of  previous  attacks 
to  insure  that  adequate  damage  has  been  inflicted.  This  survey  could  involve 
the  low-altitude  fly-over  of  ground  areas  contaminated  by  neutron  activation 
of  surface  material  and  fallout  of  radioactive  dust  and  weapon  debris  from  the 
radioactive  cloud.  This  report  addresses  this  deficiency  by  developing  an 
engineering  model  of  the  ground  accumulation  of  fallout  and  by  calculating  the 
ionizing  doses  (and  dose  rates)  which  an  aircraft  above  this  material  would 
accumulate. 

One  additional  area  is  also  investigated,  i.e.,  the  dose  rate  as  a function 
of  horizontal  distance  from  the  radioactive  area  of  ground  and  from  the  radio- 
active dust  cloud.  These  results  may  prove  useful  to  nuclear  survivability/ 
vulnerability  (S/V)  analysts  investigating  the  feasibility  of  using  radiation 
sensors  to  detect  radioactive  clouds  and  surface  areas  at  distances  sufficient 
for  avoidance  maneuvers. 
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SECTION  II 
CALCULATIONS 


SURFACE  ACTIVITY  CALCULATION 

The  calculations  of  the  ionizing  dose,  D,  accunxilated  by  an  aircraft  flying 
over  an  infinite*  radiating  plane  and  the  associated  ionizing  dose  rate,  5, 
are  conceptually  quite  simple.  The  governing  equations  can  be  derived  from 
the  geoiretries  depicted  in  figure  1. 

Ag(t)  ^ XdXde 

4tt  (h*  + X'-)  (1) 

D(t)  “ I D(t)  dt 

•'t  (2) 

where  A^  is  the  surface  activity  in  photons/cm^-hr. , h is  the  aircraft  altitude 
in  meters,  t is  time  in  hours  after  detonation,  and  \i'  is  the  atmospheric  ab- 
sorption factor,  (p''  “ 6.767  x 10"^  m“*  at  1000  meters  above  sea  level). 

Although  conceptually  simple,  the  above  calculations  in  practice  are  almost 
impossible  to  perform  because  of  the  difficulty  in  defining  the  surface  activity. 
The  above  equations  are  based  on  the  implicit  assumption  that  the  surface 
activity  is  the  result  of  uniformly  distributed  radioactive  material  and  is 
therefore  only  a function  of  time.  In  reality,  the  surface  activity  is  a com- 
plicated function  and  varies  with  location  as  well  as  time.  There  are  "hot 
spots"  and  gradients  caused  by  winds,  surface  irregularities,  detonation 
geometries,  nonuniform  cloud  fallout,  and  numerous  other  factors.  To  obtain 
even  approximate  results,  simplifying  assumptions  must  be  made.  The  major  ones 


*Infinite  in  this  context  Implies  a radius,  Rp,  such  that  the  majority  of 
photons  emitted  from  sources  outside  a sphere  of  radius  Rp,  are  absorbed 
by  the  atmosphere.  For  an  altitude  of  30,000  feet,  Rp  is  about  3000  meters, 
and  for  an  altitude  of  1000  meters,  Rp  is  about  1100  meters. 
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made  in  this  report  are  that  the  surface  activity  is  only  a function  of  time 
and  that  it  can  be  approximated  by  using  the  simple  engineering  cloud  model 
developed  in  reference  8.  This  uniform  distribution  assumption  can  be  argued 
to  yield  representative  results  for  a flyover  of  the  contaminated  area  because 
the  integrated  results  for  a flyover  should  be  similar  in  value  for  a uniform 
distribution  of  radioactive  material  as  for  the  complicated,  real  distribution. 
However,  local  results  could  vary  by  large  margins. 

The  surface  activity,  A^(t),  of  the  surface  is  assumed  to  be  a result  of 
two  mechanisms,  neutron  activation  of  surface  materials  and  fallout  of  radio- 
active material.  At  early  times,  i.e.,  less  than  10  minutes  after  detonation, 
there  is  so  much  turbulence  and  large  airborne  crater  ejecta  that  aircraft  could 
not  survive  near  the  detonation  point.  Therefore  this  analysis  shall  limit 
consideration  to  times  greater  than  10  minutes  after  detonation.  (The  same 
assumption  was  made  in  reference  7 in  the  investigation  of  radioactive  cloud 
penetrations.)  At  10  minutes  after  detonation,  the  larger  crater  ejecta  and 
other  surface  material  swept  up  into  the  cloud  during  and  immediately  after 
detonation  should  have  fallen  out  and  accumulated  on  the  ground.  Associated 
with  this  material  and  the  neutron  activetad  surface  material  is  some  activity, 
A^i  (t),  termed  the  10  minute  surface  activity.  The  remainder  of  the  radio- 
active material  is  suspended  in  the  atmosphere  and  forms  the  radioactive  dust 
cloud.  As  the  time  after  detonation  increases,  the  material  in  the  cloud  falls 
out  of  the  cloud  and  accumulates  on  the  surface.  This  cloud  iallout  material 
results  in  additional  surface  activity,  A^^(t),  termed  fallout  surface  activity. 
The  total  surface  activity  at  time  t is  the  sum  of  the  10  minute  and  fallout 
acti vi ties. 


As(t)  = Aj^(t)  + A^^(t)  (3) 

The  amount  of  radioactivity  associated  with  the  very  large  particles  (those 
which  have  impacted  the  surface  by  10  minutes  after  detonation)  and  with  neutron 
activated  surface  material  is  significant,  especially  at  early  times,  and  must 
be  considered.  However,  the  majority  of  the  radioactive  material  is  contained 
in  the  radioactive  cloud. 

For  conservation,  it  was  assumed  in  reference  a that  at  10  minutes  , aVl^  of 
the  radioactivity  resulting  from  the  detonation  was  contained  in  the  cloud  as 
coatings  on  dust  particles  or  as  condensed  weapon  debris.  In  this  analysis. 
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it  is  postulated  that  the  total  surface  radioactivity  at  10  minute-j  is  10  per- 
cent of  the  total  cloud  radioactivity  at  10  minutes.  The  only  time  dependence 
exhibited  by  the  10  minute  surface  activity  then  should  be  the  standard  one 
used  in  G'.;sstone  (ref.  14) „ 

Total  activity  in  the  cloud  is  just  the  specific  activity  of  the  cloud.  A, 
in  photops/gram  (dust)  - hour,  multiplied  by  the  total  dust  mass  in  the  cloud, 
i.e.,  A • surface  area  • cloud  height  • cloud  dust  density.  The  10  minute  sur 
face  activity  then  is  this  total  cloud  activity  multiplied  by  10  percent  f^nd 
then  divided  by  the  surface  area.  Performing  these  operations  yields  the 
relation 


^ = 1.75t"*-^ 


(4) 


where  = 100  h^a^A^  photons/hour  - centimeter^,  = 3.46  x 10"®L^.  grams 
(dust)/centimeters\  Aj  = 4.22  x 10*^  photons/gram  (dust)-hour,  is  the  load 
factor  in  megatons  of  yield  per  kilometers^,  and  h^  is  the  height  of  the  cloud 
in  meters. 

The  fallout  surface  activity,  by  definition,  is  zero  at  10  minutes  after 
detonation.  However,  as  radioactive  material  falls  out  of  the  cloud  and 
accumulates  on  the  ground,  the  fallout  surface  activity  level  begins  to  acquire 
significance.  The  fallout  surface  activity  is  a complex  function  of  time.  It 
tends  to  increase  with  time  because  radioactive  material  is  continually  accumu- 
lating on  the  surface.  But,  at  the  same  time,  it  tends  to  decrease  because  of 
the  standard  radioactive  decay  process. 

To  estimate  the  behavior  of  the  fallout  activity  an  engineering  model 
representative  of  the  after-10-minute  surface  fallout  accumulation  is  needed. 
Since  the  cloud  model  developed  in  reference  8 is  judged  to  be  reasonably 
representative  of  the  actual  cloud,  it  will  be  a basis  for  the  surface  ^allouc 
model.  The  cloud  model  is  based  on  several  assumptions.  The  major  ones  are 
100  percent  fusion  yield  of  the  warheads,  1 megaton  of  dust  lofted  per  megaton 
of  yield,  zero  wind,  fallout  of  dust  particles  based  on  Stokes  Law  type 
behavior  of  small  particles,  representative  soil  at  the  detonation  site,  and 
postulated  specific  activities  based  on  particle  size.  The  cloud  model, 
although  simple  conceptually,  was  forced  to  agree  with  the  output  of  complex 
computer  codes  (see  reference  8)  in  terms  of  dust  density  as  a function  of  time. 
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and  hence  should  be  indicative  of  the  actual  cloud.  Another  significant 
assumption  is  that  the  eaj;ly  time  turbulence  results  in  a uniform  mixture  of 
particle  sizes  over  the  pertinent  surface  area  to  an  altitude  of  15  kilometers. 
For  single  detonations,  the  pertinent  surface  area  is  a circle  with  radius  of 
kilometers,  where 

0 . 4 

R^  = 5.7W 

where  W is  the  yield  in  megatons.  For  a massive,  attack,  i.e.,  many  detonations 
over  a relatively  small  area,  the  cloud  radii  may  overlap,  and  the  total  sur- 
face area  must  be  estimated  by  other  techniques. 

The  cloud  is  characterized  by  two  parameters:  the  dust  density,  p^,  in 

grams  (dust)  per  centimeter^  and  the  specific  activity.  A,  in  photons  per 
gram  (dust)-hour.  These  relations  are  expressed  in  the  form 


ij(r,t)  - 3.42  X 10‘^  r-®-'^ 

O.lp  < 

r < R'(t) 

(5) 

A(r,t)  = 5 X 10“ 

O.lp  < 

r < 20u 

A(r,t)  « 10“ 

20u  < 

r < 10-p 

(6) 

where  R^(t)  = particle  size  fallout  function  in  microns 

= 5.04  (t"*‘®  + 1.313  t“°*’)  + 0.316  , r is  the  particle  radius 

L J 

in  microns,  and  t is  the  time  after  detonation  in  hours. 

For  the  fallout  activity  model,  the  10  minute  cloud  will  be  assumed  to  be 
uniformly  distributed  over  the  volume  discussed  previously.  All  of  the  parti- 
cles in  the  cloud  are  falling  with  some  vertical  velocity,  V.  It  is  assumed 
the  particles  are  spherical  and  that  there  are  two  vertical  forces  acting  on 
each  particle:  the  weight  of  the  particle  and  the  aerodynamic  drag.  Figure  2 

shows  these  forces.  Summing  vertical  forces  yields  the  relationship 

^ Ttr^  g = J-  r CjjTTr^ 

Solving  for  the  vertical  velocity  yields 
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Ft=  1/2  Po  Cd 


DRAG  FORCE 


W : 


f 


Figure  2.  Vertical  Forces  Acting  on  a Oust  Particle  in  the  Dust  Cloud 

where  g is  the  gravitational  constant,  is  the  density  of  the  particle, 
is  the  density  of  the  air,  r is  the  particle  radius,  and  Cp  is  the  drag  coef- 
ficient of  the  sphere.  This  drag  coefficient  has  been  experimentally  measured, 
and  its  behavior  as  a function  of  Reynolds  Number*  is  depicted  in  figure  3, 
which  is  adapted  from  Binder  (ref.  12).  This  curve  was  approximated  over  the 
various  Reynolds  Number  ranges  depicted  in  the  figure  by  the  following  expres- 
sions. (The  approximations  are  the  dashed  lines.) 

Cp  = 24/Re  O.lp  < r < 60p  (8) 


*The  Reynolds  Number,  Re,  is  a nondimensional  parameter  and  is 


where  d is  the  diameter  of  the  sphere  and  p.  the  dynamic  viscosity  of  air. 

a 
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Cq  = 14. 79/Re'’- ^ 60u  < r < IZOop 


= 0.4557 


1200u  < r < 10^1 


These  expressions  in  conjunction  with  the  previous  equation  define  the  particle 
settling  velocity  as  a function  of  particle  size.  Note  that  large  particles 
have  larger  velocities  and  hence  accumulate  faster  on  the  surface. 

Now  consider  a portion  of  the  cloud  over  a square  centimeter  of  the  surface. 
The  particles  in  this  volume  at  various  times  are  depicted  in  particle  size- 
altitude  space  in  figure  4.  At  10  minutes,  particles  are  uniformly  distributed 
over'  all  altitudes.  Since  the  larger  particles  fall  faster,  the  higher  altitude 

t s 10  MINUTES 


♦3  > *2  > ♦•  > »0  MINUTES 


10  10' 
PARTICLE  RADIUS 


Figure  4.  Particle  Size  Distribution  with  Altitude  and  Time 
for  Radioactive  Dust  Clouds 
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dust  cloud  strata  will  be  depleted  of  these  particles  at  late  times.  The  total 
mass  of  particles  in  the  depleted  regions  can  be  calculated  because  the  vertical 
velocity  of  each  particle  size  is  known.  The  dust  mass  multiplied  by  the 
specific  activity  (both  are  functions  of  particle  size)  and  then  integrated  over 
the  depleted  region  in  the  depleted  space  indicated  in  figure  4,  results  in 
the  surface  fallout  activity, 

A^^(t)  = j I A(r,t)  dhdr 

Jo.  UiJh3(r,t) 

1200m  h 

+ I I <’d^’"’''o^  A(r,t)  dhdr 

JsOp  Jh^(r,t) 

iQ.OOOu  hg 

+ 1 I ‘^d^**’^o^  A(r,t)  dhdr 
,/i200m  *'hj(r,t)  (11) 

where  t^  is  10  minutes  after  detonation  and  is  constant,  h^  is  the  cloud  height 
(assumed  to  be  15,000  meters),  and  from  the  particle  velocity  relations 

h,  (r,t)  «=  h^  - 1797.6  r^'^t-tg)  1200u  ••  r ■ I0"ii 

h_,  (r,t)  = - 51.9  r (t-t^)  60u  < ‘ * 1200m  ' 

h^  (r,t)  * h^  - 0.86  r^ 

Performing  the  indicated  operations  yields  the  analytical  expression 

= 1.0694  1.787  f-'  (^2) 

The  integrations  were  performed  on  the  Air  Force  Weapons  Laboratory  CDC  6600 
computer.  The  results  are  depicted  in  figures  5 and  6.  Figure  5 shows  the 
generalized  results  in  terms  of  Kg,  and  figure  6 shows  the  specific  results  for 
a 4-megaton  surface  detonation.  The  results  for  a specific  threat  provide 
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physical  insight  into  the  levels  of  activity  that  may  be  observed  from  a typical 
detonation.  Note  that  for  times  less  than  1 hour  after  detonation  the  IJ  min- 
ute surface  activity  is  the  dominant  factor.  However,  at  late  times  it  becomes 
less  significant,  and  the  fallout  surface  activity  becomes  dominant  (figs . 5 & 6). 

Note  that  in  case  the  cloud  is  relocated  by  the  wind,  only  the  fallout  sur- 
face activity  would  be  used  to  calculate  doses  and  dose  rates  for  flyovers  of 
the  new  surface  where  A^j(t)  - 0.  The  dose  rate  and  dose  at  an  aircraft  at 
altitude  h would  be  determined  using  equations  (1)  and  (2)  with  A (t)  = A (t). 

FLYOVER  DOSE  AND  DOSE  RATE  CALCULATIONS 

The  calculations  of  the  ionizing  dose  accumulated  by  an  aircraft  flying 
over  the  contaminated  surface  at  an  altitude  h and  of  the  ionizing  dose  rate  as 
a function  of  time  after  detonation  can  now  be  made.  The  governing  equations 
are  equations  (1)  and  (2),  introduced  in  the  previous  section,  with  the  surface 
activity,  A^,  in  these  equations  being  defined  by  equations  (3),  (4),  and  (11). 
Performing  the  indicated  operations  for  a baseline  aircraft  altitude  of  100 
meters  above  the  contaminated  surface  (assumed  to  be  1000  meters  above  sea 
level)  yields  a dose  rate,  D^(t)  in  aircraft  defined 

by  the  equation 


D,(t) 

= 2.1  t*''-''’  -6.7  X 10"’  f'-' 

where  = C and  C = 4.88  x 10" rads  (tissue)/photon/cm^ 

(for  the  assumed  1 MeV  photons). 

The  dose  0|(t)  in  rads  (tissue)  accumulated  by  the  crew  of  the  aircraft  from 
some  contaminated  area  entry  time  t^  and  to  some  exit  time  t (where  both 
times  are  in  units  of  hours  after  detonation)  is  defined  by  the  equation 

Df(t) 

— = 2.61  (t*’-"  - t^“-®)  + 3.36  X lO"-"  (t'“*=^  - t^“*’-‘’)  (14) 

The  dose  rate  results  are  presented  In  figures  7 to  10  and  the  dose  results 
in  figures  11  to  14.  The  generalized  results  in  figures  7,  8,  11,  and  12  are 
presented  in  terms  of  the  constant,  K,  rads  (tissue)/hour.  This  constant  is 
used  to  generalize  the  results.  For  a given  threat  this  constant  is  fixed. 
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and  specific  dose  rates  and  doses  can  be  obtained.  Figures  9,  10,  13,  and  14 
correspond  to  a specific  threat  and  are  presented  to  give  the  user  a "feel"  for 
the  type  of  results  he  may  encounter.  These  figures  are  based  on  a 4 megaton 
(MT)  surface  detonation,  i.e.,  = 0.0129  MT/km^. 

All  of  these  results  are  for  a baseline  aircraft  altitude  of  100  meters. 
Scaling  factors  to  make  these  results  applicable  to  other  altitudes  are  pre- 
sented in  figure  15.  To  use  this  scaling  factor,  first  determine  the  altitude 
of  interest,  then  determine  the  pertinent  scaling  factor  from  figure  15.  (For 
example,  if  the  aircraft  altitude  is  200  meters,  the  scaling  factor  is  0.34.) 

The  scaling  factor  is  applied  to  the  ordinates  or  figures  7 to  14  by  simple 
multiplication. 

Note  that  the  dose  results  are  presented  as  families  of  curves,  one  curve 
for  each  contaminated  area  entry  time.  Therefore,  for  each  flyover,  the  threat, 
the  aircraft  altitude,  the  entry  time,  and  the  contaminated  area  exit  time 
must  be  known  in  order  to  "fix"  the  ionizing  dose  accumulated  by  the  aircraft. 
The  instantaneous  ionizing  do*,e  rate  for  any  time  during  flyover  can  also  be 
"fixed."  Note  that  the  dose  and  dose  rate  results  are  presented  in  rads 
(tissue)  which  are  directly  applicable  to  the  aircrew.  If  the  results  are 
desired  in  terms  of  ruds  (silicon)  for  application  to  electronic  equipment,  then 
the  conversion  factor  to  use  (for  the  assumed  1 MeV  photons)  is 

rads  (silicon)  0.922  rads  (tissue)  (15) 

It  is  noted  that  during  the  flyover  of  the  contaminated  area,  the  aircraft 
may  also  be  penetrating  the  radioactive  cloud  above  this  area.  Therefore,  the 
total  ionizing  dose  accumulated  would  be  the  sum  of  the  doses  from  (1)  the 
contaminated  surface  and  (2)  the  radioactive  cloud.  The  first  dose  can  be 
estimated  using  the  methods  and  results  presented  in  this  paper.  The  second 
dose  can  be  estimated  by  using  the  techniques  and  methods  presented  in 
reference  8. 

To  illustrate  the  technique  of  using  these  results,  an  example  will  be 
presented.  Assume  that  a 4 megaton  weapon  had  detonated  on  the  surface 

R = 5.7(4)®*“ 

R = 9.92  kilometers 


25 


AFWL-TR-75-214 


and 

, , jl.me3atgns „ HJl 

' n(9.92)’  km'  km' 

Assume  that  an  aircraft  is  flying  over  the  radioactive  surface  at  an  altitude 
of  150  meters  at  a velocity  of  200  meters/second.  At  2 hours  after  detonation 
the  aircraft  enters  the  contaminated  area  to  conduct  a damage  assessment  sur- 
vey. Since  its  velocity  is  200  meters/second  and  the  diameter  of  the  contami- 
nated area  is  about  2 x lO**  meters,  the  aircraft  exits  the  radioactive  area 
1.7  minutes  after  entry. 

With  a contaminated  area  entry  time,  t. , and  exit  time,  t^,  the  generalized 
dose  results  of  figures  11  and  12  could  be  used.  And  for  relatively  long 
flyovers,  i.e.,  (t^  - t^)  > 10  minutes,  reasonably  accurate  results  could  be 
obtained.  However,  for  short  flyovers,  tlie  i jsolution  is  poor,  and  the  results 
obtained  from  the  dose  curves  are  questionable.  Better  estimates  for  short 
flyovers  can  be  approximated  by  using  the  dose  rate  results  of  figures  9 and 
10.  The  dose  is  simply 

D ~ D (t^  - t^)  (16) 

where  D is  the  dose  rate  at  the  time  (t^  + t.j  )/2  and  the  times  are  in  hours. 
Using  this  technique,  the  flyover  dose  accumulated  by  the  aircraft  for  the 
situation  in  the  previous  paragraph  is  (244  xO.55)  rads(  ti ssue)/hour  (1.7/60) 
hours,  or  3.8  rads(tissue) , where  0.55  is  the  scaling  factor  obtained  from 
figure  15,  and  244  is  the  dose  rate  at  2 hours  after  detonation  obtained  from 
figure  9. 

IONIZING  DOSE  RATE  AT  VARIOUS  HORIZONTAL  DISTANCES  FTOM  RADIOACTIVE  SURFACE 
AREAS  AND  RADIOACTIVE  DUST  CLOUDS 

An  aircraft  engaging  in  conducting  damage  assessment  surveys  of  targets 
previously  attacked  may  be  required  to  fly  over  radioactive  surface  areas  to 
accomplish  its  mission.  However,  to  keep  the  ionizing  dose  accumulation  to 
acceptable  levels,  avoidance  of  radioactive  clouds  and  radioactive  surface 
areas  not  specifically  designated  as  its  responsibility  may  be  desirable.  A 
prerequisite  to  avoidance  is  detection.  An  important  consideration  in  the 
design  and  development  of  hardware  to  accomplish  detection  of  radioactive  dust 
clouds  and/or  surface  areas  is  the  determination  of  ionizing  dose  rates  as  a 
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function  of  horizontal  distance  from  the  boundary  of  these  radioactive  clouds/ 
# Fjrface  areas.  Sensors  should  be  developed  sensitive  enough  to  allow  detection 

and  avoidance  maneuvers  prior  to  penetration  of  the  cloud  and/or  flyover  of 
the  surface  area  if  this  technique  is  to  be  of  optimum  utility. 

i 

! Consider  first  the  radioactive  cloud.  It  is  assumed  that  the  cloud  boundary 

is  well  defined,  that  all  the  airborne  radioactive  material  is  uniformly  dis- 
tributed throughout  the  cloud,  and  that  the  cloud  boundary  is  a vertical  plane 
perpendicular  to  tne  aircraft  flightpith.  This  situation  is  schematically 
depicted  in  figure  16.  This  figure  shows  the  aircraft  at  some  distance,  d, 
from  the  cloud  boundary  and  the  ring  volume  element  of  the  cloud.  The  ionizing 
dose  rate  at  the  aircraft  due  to  the  radioactive  dust  particles  in  this 
j element  is 

I r(t) 

D^(t)  = Pjj(r,t)  A{r,t)  dr  (17) 

Jo. Ip 

The  dose  rate  at  the  aircraft  due  to  all  of  the  particles  in  the  cloud  then  is 
-n/2  --  -R'(t) 

D^(t)  = I I I P^(r.t)  A(r,t)  Ce“^  ^ sin  0 drdRdO 
Jo  JdsecG  Jo.lp 

The  integration  was  performed  numerically,  and  the  dose  rates  as  a function 
of  the  time  after  detonation  are  shown  in  figures  17  and  18  for  a baseline 
distance  of  4000  meters  between  the  aircraft  and  the  cloud  boundary.  Tigure 
19  provides  scaling  factors  if  the  dose  rate  corresponding  to  another  distance 
is  required.  Note  that  this  scaling  factor  is  1.0  for  a distance  of  4000 
meters.  If  the  distance  is  less,  the  scaling  factor  is  larger  and  vice  versa. 
This  scaling  factor  is  used  to  scale  the  ordinate  of  the  dose  rate  graph  of 
figures  17  and  18. 

The  radioactive  surface  area  situation  is  depicted  in  figure  20.  The  air- 
craft is  flying  at  an  altitude  of  h meters  and  is  a distance  d meters  from 
the  boundary  of  the  radioactive  surface  area.  The  dose  rate  at  the  aircraft 

I from  the  radioactive  material  in  the  surface  area  element  is 
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Figure  19.  Scaling  Factor  for  Distance  from  Cloud  Boundary 
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Figure  20.  Aircraft  Approaching  an  Infinite  Radioactive  Surface  Area 
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Oj(t) 


CA.(t)  ^ 


4I1(R^^  + h^) 


dRj  dG 
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The  dose  rate  at  the  aircraft  due  to  all  of  the  radioactive  material  on  the 
surface  is 


nn 

“sf'l*  1 

Jo  Jd 


CAM)  e'^'  ^^^s' 


ds  sec  0 


411  (Rj.2  + hO 


dRjdG 


(20) 


Performing  the  indicated  operations,  the  following  expression  is  obtained. 

DOt)  9.76  X 10“‘  A.(t) 

- , , „ e , , , , , ■ , — I ,,  — 

^ (21) 


For  a baseline  altitude  of  100  meters  and  a horizontal  distance  of  1000 
meters,  the  dose  rate  is  shown  in  figures  21  to  24  as  a function  of  time  after 
detonation.  For  other  altitudes  and  horizontal  distances,  scaling  factors 
from  figure  25  must  be  used  to  scale  the  ordinate  of  figures  21  to  24. 

A specific  aircraft  being  analyzed  has  well  defined  performance  character- 
istics. Therefore,  a horizontal  distance  can  be  defined  which  corresponds  to 
successful  avoidance  of  the  cloud  and/or  surface  area.  For  this  distance, 
an  ionizing  dose  rate  can  be  defined  from  the  results  above.  An  onboard  radia- 
tion sensor  must  be  capable  of  detecting  an  instantaneous  dose  rate  of  that 
magnitude  for  the  avoidance  technique  to  be  feasible. 
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